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Abstract. In this paper we present applications of methods from wavelet analysis to 

polynomial approximations for a number of accelerator physics problems. According 

to variational approach in the general case we have the solution as a multiresolution 

(multiscales) expansion in the base of compactly supported wavelet basis. We give 

Oh| extension of our results to the cases of periodic orbital particle motion and arbitrary 

variable coefficients. Then we consider more flexible variational method which is based 

^vq | on biorthogonal wavelet approach. Also we consider different variational approach, 

£> ■ which is applied to each scale. 
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This is the first part of our two presentation in which we consider applications of 
methods from wavelet analysis to nonlinear accelerator physics problems. This is a 
continuation of results from [l]-[6], which is based on our approach to investigation 
of nonlinear problems - general, with additional structures (Hamiltonian, symplec- 
tic or quasicomplex) , chaotic, quasiclassical, quantum, which are considered in the 
framework of local (nonlinear) Fourier analysis, or wavelet analysis. Wavelet anal- 
ysis is a relatively novel set of mathematical methods, which gives us a possibility 
to work with well-localized bases in functional spaces and with the general type of 
operators (differential, integral, pseudodifferential) in such bases. 

We consider application of multiresolution representation to general nonlinear dy- 
namical system with polynomial type of nonlinearities. In part II we consider this 
very useful approximation in the cases of orbital motion in storage rings, a particle 
in the multipolar field, effects of insertion devices on beam dynamics, spin orbital 
motion. Starting in part III A from variational formulation of initial dynamical 
problem we construct via multiresolution analysis (part III B) explicit representa- 
tion for all dynamical variables in the base of compactly supported (Daubechies) 
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wavelets. Our solutions (part III C) are parametrized by solutions of a number 
of reduced algebraical problems one from which is nonlinear with the same degree 
of nonlinearity and the rest are the linear problems which correspond to particu- 
lar method of calculation of scalar products of functions from wavelet bases and 
their derivatives. Then we consider further extension of our previous results. In 
part V we consider modification of our construction to the periodic case, in part 
VI we consider generalization of our approach to variational formulation in the 
biorthogonal bases of compactly supported wavelets and in part VII to the case 
of variable coefficients. In part IV we consider the different variational approach 
which is based on ideas of para-products (A) and approximation for multiresolution 
approach, which gives us possibility for computations in each scale separately (B). 

II PROBLEMS AND APPROXIMATIONS 

Below we consider a number of examples of nonlinear accelerator physics prob- 
lems which are from the formal mathematical point of view not more than nonlinear 
differential equations with polynomial nonlinearities and variable coefficients. 

A Orbital Motion in Storage Rings 

We consider as the main example the particle motion in storage rings in standard 
approach, which is based on consideration in [7]. Starting from Hamiltonian, which 
described classical dynamics in storage rings 

H{f,P,t) = c{7i 2 + m 2 c 2 } 1 / 2 + e<f) (1) 

and using Serret-Frenet parametrization, we have the following Hamiltonian for 
orbital motion in machine coordinates: 

H{x,p x ,z,p z ,a,p a ;s) =p a -[l + f(p ff ] ■ [1 + K x ■ x + K z ■ z] x (2) 
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Then, after standard manipulations with truncation of power series expansion 
of square root we arrive to the following approximated Hamiltonian for particle 
motion: 
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and the corresponding equations of motion: 
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Then we use series expansion of function f(p a ) from [7]: 
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and the corresponding expansion of RHS of equations (4). In the following we 
take into account only an arbitrary polynomial (in terms of dynamical variables) 
expressions and neglecting all nonpolynomial types of expressions, i.e. we consider 
such approximations of RHS, which are not more than polynomial functions in 
dynamical variables and arbitrary functions of independent variable s ("time" in 
our case, if we consider our system of equations as dynamical problem). 



B Particle in the Multipolar Field 

The magnetic vector potential of a magnet with In poles in Cartesian coordinates 



is 



A = J2Knfn(x,y): 



(6) 



where /„ is a homogeneous function of x and y of order n. 
The real and imaginary parts of binomial expansion of 

fn(x,y) = (x + iy) n 



(7) 



correspond to regular and skew multipoles. The cases n = 2 to n = 5 correspond 
to low-order multipoles: quadrupole, sextupole, octupole, decapole. 

Then we have in particular case the following equations of motion for single par- 
ticle in a circular magnetic lattice in the transverse plane (x, y) ([8] for designation): 
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and the corresponding Hamiltonian: 
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Then we may take into account arbitrary but finite number in expansion of RHS of 
Hamiltonian (9) and from our point of view the corresponding Hamiltonian equa- 
tions of motions are also not more than nonlinear ordinary differential equations 
with polynomial nonlinearities and variable coefficients. 



C Effects of Insertion Devices on Beam Dynamics 

Assuming a sinusoidal field variation, we may consider according to [9] the ana- 
lytical treatment of the effects of insertion devices on beam dynamics. One of the 
major detrimental aspects of the installation of insertion devices is the resulting 
reduction of dynamic aperture. Introduction of non-linearities leads to enhance- 
ment of the amplitude-dependent tune shifts and distortion of phase space. The 
nonlinear fields will produce significant effects at large betatron amplitudes. 



The components of the insertion device magnetic field used for the derivation of 
equations of motion are as follows: 

k 
B x = -j- • B sinh(k x x) smh(k y y) cos(kz) 

ky 

By = B cosh(k x x) cosh(k y y) cos(kz) (10) 
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B z = ——B cosh(k x x) sinh(k y y) sm(kz) 
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with k 2 + k 2 = k 2 = (27r/A) 2 , where A is the period length of the insertion device, -Bo 
is its magnetic field, p is the radius of the curvature in the field Bq. After a canonical 
transformation to change to betatron variables, the Hamiltonian is averaged over 
the period of the insertion device and hyperbolic functions are expanded to the 
fourth order in x and y (or arbitrary order). 
Then we have the following Hamiltonian: 
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We have in this case also nonlinear (polynomial with degree 3) dynamical system 
with variable (periodic) coefficients. As related cases we may consider wiggler and 
undulator magnets. We have in horizontal x — s plane the following equations 

x = -s—B z (s) (12) 

7727 

e 
s = x B z (s), 

where magnetic field has periodic dependence on s and hyperbolic on z. 

D Spin-Orbital Motion 

Let us consider the system of equations for orbital motion 

dq _ dH orb dp __ dH orb 

dt dp dt dq 

and Thomas-BMT equation for classical spin vector (see [10] for designation) 
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? — (^i? 92, 93), P = (pi,P2,P3) are canonical position and momentum, s = (si, s 2 , s 3 ) 
is the classical spin vector of length K/2, n = (tti, 7r 2 , ^3) is kinetic momentum 
vector. We may introduce in 9-dimensional phase space z = (q, p, s) the Poisson 
brackets 

{f(z),g(z)} = f q g p - f p g q + [f s x g s \ ■ s (16) 

and the corresponding Hamiltonian equations: 

ch 

d~t ={Z > Hh (17) 

with Hamiltonian 

H = H orb (q,p,t)+w(q,p,t)- s. (18) 

More explicitly we have 

dq dH orb d(w ■ s) 



dt dp dp 

dp _ dH orb _ d(w ■ s) 

dt dq dq 

ds 

— = \w X s\ 

dt L J 

We will consider this dynamical system also in our second paper in this volume via 

invariant approach, based on consideration of Lie-Poison structures on semidirect 

products of groups. 

But from the point of view which we used in this paper we may consider the 

similar approximations as in preceding examples and then we also arrive to some 

type of polynomial dynamics. 

Ill POLYNOMIAL DYNAMICS 

The first main part of our consideration is some variational approach to this 
problem, which reduces initial problem to the problem of solution of functional 
equations at the first stage and some algebraical problems at the second stage. We 
have the solution in a compactly supported wavelet basis. Multiresolution expan- 
sion is the second main part of our construction. The solution is parameterized 
by solutions of two reduced algebraical problems, one is nonlinear and the second 
is some linear problem, which is obtained from one of the next wavelet construc- 
tions: Fast Wavelet Transform (FWT), Stationary Subdivision Schemes (SSS), the 
method of Connection Coefficients (CC). 



A Variational Method 

Our problems may be formulated as the systems of ordinary differential equations 

fci/dt = fi( Xj ,t), (i,j = l,...,n) (20) 

with fixed initial conditions £j(0), where fa are not more than polynomial functions 
of dynamical variables Xj and have arbitrary dependence of time. Because of time 
dilation we can consider only next time interval: < t < 1. Let us consider a set 
of functions 

$i(t) = xdji/dt + f iVi (21) 

and a set of functionals 

Fi(x) = I $i(t)dt - x iVi |S, (22) 



where yi(t)(i/i(0) = 0) are dual variables. It is obvious that the initial system and 
the system 

F(x) = (23) 

are equivalent. In part VI we consider more general approach, which is based on 
possibility taking into account underlying symplectic structure and on more useful 
and flexible analytical approach, related to bilinear structure of initial functional. 
Now we consider formal expansions for Xi,yf. 

Xl (t)= Xl (0) + J2^Mt) %(t) = ]>»(*), (24) 

k r 

where because of initial conditions we need only y?fc(0) = 0. Then we have the 
following reduced algebraical system of equations on the set of unknown coefficients 
A^ of expansions (24): 

$>*rA? - 7i(A,-) = (25) 

k 

Its coefficients are 

Vkr= <p' k (t)(p r (t)dt, Yi = fi(Xj,t)(p r (t)dt. (26) 

JO Jo 

Now, when we solve system (25) and determine unknown coefficients from formal 
expansion (24) we therefore obtain the solution of our initial problem. It should be 
noted if we consider only truncated expansion (24) with N terms then we have from 
(25) the system of N x n algebraical equations and the degree of this algebraical 
system coincides with degree of initial differential system. So, we have the solution 
of the initial nonlinear (polynomial) problem in the form 



N 

x i (t)=x i (0) + ^2\!X k (t), (27) 
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where coefficients Xf are roots of the corresponding reduced algebraical problem 
(25). Consequently, we have a parametrization of solution of initial problem by 
solution of reduced algebraical problem (25). The first main problem is a problem 
of computations of coefficients of reduced algebraical system. As we will see, these 
problems may be explicitly solved in wavelet approach. 

Next we consider the construction of explicit time solution for our problem. The 
obtained solutions are given in the form (27), where X^{t) are basis functions and 
\ l k are roots of reduced system of equations. In our first wavelet case Xj,(t) are 
obtained via multiresolution expansions and represented by compactly supported 
wavelets and A^ are the roots of corresponding general polynomial system (25) with 
coefficients, which are given by FWT, SSS or CC constructions. According to the 
variational method to give the reduction from differential to algebraical system of 
equations we need compute the objects 7^ and fiji, which are constructed from 
objects: 

(Ji= Xi(r)dr, 
Jo 

Vij = / X i (r)X j (r)dr, 

Hji= /x;(r)X,(r)dr, (28) 



Pkij= I'' X^X^X^dr 

J u 



for the simplest case of Riccati systems, where degree of nonlinearity equals to two. 
For the general case of arbitrary n we have analogous to (28) iterated integrals 
with the degree of monomials in integrand which is one more bigger than degree of 
initial system. 

B Wavelet Framework 

Our constructions are based on multiresolution approach. Because affine group 
of translation and dilations is inside the approach, this method resembles the action 
of a microscope. We have contribution to final result from each scale of resolution 
from the whole infinite scale of spaces. More exactly, the closed subspace Vj(j G Z) 
corresponds to level j of resolution, or to scale j. We consider a r- regular mul- 
tiresolution analysis of L 2 (R n ) (of course, we may consider any different functional 
space) which is a sequence of increasing closed subspaces Vf 

...V-2 C V.! C V C V 1 C V 2 C ... (29) 

satisfying the following properties: 



n^- = o, \Jv j = l 2 (bp), 

jez jez 

f(x) e Vj <=> f(2x) e V j+1 , 

f(x)ev <=>f(x-k)ev , ykez n . (30) 

There exists a function ip G V such that {ip 0tk (x) = ip(x — k),k G Z n } forms a 
Riesz basis for V . 

The function p is regular and localized: (p is C r_1 , c^^ 1 ) is almost everywhere 
differentiable and for almost every x G R n , for every integer a < r and for all 
integer p there exists constant C p such that 

\d a <p(x)\<C p (l + \x\)- p (31) 

Let (f(x) be a scaling function, ip(x) is a wavelet function and </?j(x) = •pix — i). 
Scaling relations that define tp, ip are 

7V-1 N-l 

fix) = £ a k¥( 2x ~ k) = J2 a k ip k {2x), (32) 

fc=0 fc=0 

N-2 

i>{x) = E (-1) Vn¥>(2* + k). (33) 

fc=-i 

Let indices £,j represent translation and scaling, respectively and 

p jt {x) = 2 j/2 p{2 j x - £) (34) 

then the set {pj tk },k G Z n forms a Riesz basis for Vj. The wavelet function ^ is 
used to encode the details between two successive levels of approximation. Let Wj 
be the orthonormal complement of Vj with respect to Vj + \\ 

V j+ i = Vj®Wj. (35) 

Then just as Vj is spanned by dilation and translations of the scaling function, so 
are Wj spanned by translations and dilation of the mother wavelet ipj k (x), where 

ip jk (x) = 2 j/2 ijj(2 j x-k). (36) 

All expansions which we used are based on the following properties: 

{^jk}, j,k G Z is a Hilbertian basis of L 2 (R) 
{</2jfc}j>o,fcez is an orthonormal basis forL 2 (R), 

oo 

L 2 (R) = V ®W 3 , (37) 

3=0 

or 

{(po,k, i>j,k}j>o,kez is an orthonormal basis forL 2 (R). 



C Wavelet Computations 

Now we give construction for computations of objects (28) in the wavelet case. 
We use compactly supported wavelet basis: orthonormal basis for functions in 
L 2 (R). 

Let be / : R — ► C and the wavelet expansion is 

oo 

f( x ) = J2 c Mt{x) + J2 J2 °jkipjk(x) (38) 

eez j=okez 

If in formulae (38) Cjk = for j > J, then f(x) has an alternative expansion 
in terms of dilated scaling functions only f(x) = J2 Cje<fije(x). This is a finite 

wavelet expansion, it can be written solely in terms of translated scaling functions. 
Also we have the shortest possible support: scaling function DN (where N is even 
integer) will have support [0, N — 1] and N/2 vanishing moments. There exists 
A > such that DN has AiV continuous derivatives; for small N, X > 0.55. To 
solve our second associated linear problem we need to evaluate derivatives of f(x) 
in terms of <p(x). Let be ip™ = d n </^(x)/dx n . We consider computation of the 
wavelet - Galerkin integrals. Let f d (x) be d-derivative of function f(x), then we 
have f d (x) = Y^e c iPe( x ) > and values tp d {x) can be expanded in terms of <p(x) 

<^(x)=]TA m y9 m (x), (39) 

m 

oo 

Ki= ip d {x)ip m (x)dx, 



-oo 



where A m are wavelet-Galerkin integrals. The coefficients A m are 2-term connection 
coefficients. In general we need to find (di > 0) 

00 

Ktt= j\\vfM)dx (40) 

— oo 

For Riccati case we need to evaluate two and three connection coefficients 

oo 
/*oo P 

Af d * = / <p d ^x)cp d e 2 (x)dx, A dM = / <p dl (x)<p d2 (x)<pt(x)dx (41) 



According to CC method [11] we use the next construction. When iV in scaling 
equation is a finite even positive integer the function tp{x) has compact support 
contained in [0, N — 1]. For a fixed triple (di, di-, d-$) only some A^ 2 3 are nonzero: 
2-N<i<N-2, 2-N<m<N-2, \£ - m\ < N - 2. There are 
M = 3N 2 — 9N+7 such pairs (£, m). Let A dld2d ' A be an M-vector, whose components 
are numbers A^ 2 3 . Then we have the first reduced algebraical system : A satisfy 
the system of equations (d = di + d 2 + d 3 ) 



AA — 2 A 1 2 3 , Ai >m - q>r — / t QpQq-2£+p a r-2m+p (42) 

P 

By moment equations we have created a system of M + d + 1 equations in M 
unknowns. It has rank M and we can obtain unique solution by combination of LU 
decomposition and QR algorithm. The second reduced algebraical system gives us 
the 2-term connection coefficients: 

AA d ld2 = 2 i-d A d ld2) d = di + ^ A ^ = £ apaq _ 2e+p (43) 

v 

For nonquadratic case we have analogously additional linear problems for objects 
(40). Solving these linear problems we obtain the coefficients of nonlinear alge- 
braical system (25) and after that we obtain the coefficients of wavelet expansion 
(27). As a result we obtained the explicit time solution of our problem in the base 
of compactly supported wavelets. We use for modelling D6, D8, D10 functions and 
programs RADAU and DOPRI for testing. 

In the following we consider extension of this approach to the case of periodic 
boundary conditions, the case of presence of arbitrary variable coefficients and more 
flexible biorthogonal wavelet approach. 

IV EVALUATION OF NONLINEARITIES SCALE BY 

SCALE 

A Para-product and Decoupling between Scales 

But before we consider two different schemes of modification of our variational 
approach in the case when we consider different scales separately. For this reason 
we need to compute errors of approximations. The main problems come of course 
from nonlinear terms. We follow approach from [12]. 

Let Pj be projection operators on the subspaces Vj,j G Z: 

Pj : L 2 (R) -> V 3 (44) 

(Pjf)(x) = Y,< f><Pj,k > VjA x ) 
k 

and Qj are projection operators on the subspaces Wf 

Q 3 = Pj^ - Pj (45) 

So, for u G £ 2 (R) we have Uj = PjU and Uj G Vj, where {Vj},j G Z is a 
multiresolution analysis of L 2 (R). It is obviously that we can represent Mq in the 
following form: 

ul = 2j2(P jU )(Q jU ) + f2(Q jU )(Q jU ) + u 2 n (46) 

3=1 3=1 



In this formula there is no interaction between different scales. We may consider 
each term of (46) as a bilinear mappings: 

M 3 VW : Vj x Wj -► L 2 (R) = V^^Wy (47) 



M 3 WW : W,- x Wj -► L 2 (R) = V$ ®j>>j Wj, (48) 

For numerical purposes we need formula (46) with a finite number of scales, but 
when we consider limits j — ► oo we have 



'U 
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E(2P i ti + Q J -«)(Q i «), (49) 
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which is para-product of Bony, Coifman and Meyer. 

Now we need to expand (46) into the wavelet bases. To expand each term in (46) 
into wavelet basis, we need to consider the integrals of the products of the basis 
functions, e.g. 



Mtf ww (k,k',£) = / ^{{xW^xW, (x)dx, (50) 

J — oo 

where j' > j and 

^{(x) = 2- j/2 *fj{2- j x-k) (51) 

are the basis functions. If we consider compactly supported wavelets then 

M^ ww (k,k'J) = for \k-k'\>k , (52) 

where ko depends on the overlap of the supports of the basis functions and 

W r www {k - k', Tk-i)\<C- 2- rXM (53) 

Let us define jo as the distance between scales such that for a given e all the 
coefficients in (53) with labels r = j — j', r > jo have absolute values less than e. 
For the purposes of computing with accuracy e we replace the mappings in (47), 
(48) by 

M 3 VW : Vj x Wj - Vj ®j< f < j0 W f (54) 

M^ w : Wj x Wj -> Vj ®j<j'< j0 Wj, (55) 

Since 

Vj ®j<f<j Wy = Vj,.-, (56) 



and 

Vj C V j0 ^, Wj C V j0 ^ (57) 

we may consider bilinear mappings (54), (55) on Vj _i x V} _i. For the evaluation of 
(54), (55) as mappings V} _i x Vj _i — ► V} _i we need significantly fewer coefficients 
than for mappings (54), (55). It is enough to consider only coefficients 

M(k,k',e) = 2- J/2 (p(x - k)<p(x - k')(p(x - £)dx, (58) 

Joo 

where tp(x) is scale function. Also we have 

M(k, k', £) = 2~ j ' 2 M (k - £, k' - £), (59) 

where 

Mo{p,q)= \ ip(x — p)<f(x — q)<p(x)dx (60) 



Now as in section (3C) we may derive and solve a system of linear equations to find 
M (p,q). 

B Non-regular Approximation 

We use wavelet function ip(x), which has k vanishing moments / x h ip(x)dx = 0, 
or equivalently x h = J2 ceipe(x) for each k, < k < K. 

Let Pj again be orthogonal projector on space Vj. By tree algorithm we have 
for any u G £ 2 (R) and £ G Z, that the wavelet coefficients of Pe(u), i.e. the set 
{< u,ipj y k >, j < £ — 1, k G Z} can be compute using hierarchic algorithms from 
the set of scaling coefficients in Vi, i.e. the set {< u,ip^k >,/c G Z} [13]. Because 
for scaling function (p we have in general only J (p(x)dx = 1, therefore we have for 
any function u G L 2 (R): 

lim | 2 j/2 < u, ip jyk > -u(x) |= (61) 

j— >oo,fc2~3— +a; 

If the integer n(y?) is the largest one such that 

I x a ip(x)dx = for l<a<n (62) 

then if u G C^ n+1 ^ with u( n+1 ^ bounded we have for j — * oo uniformly in k: 

| 2 j/2 < u, (p jtk > -u{k2~ j ) |= 0(2" j(n+1) ). (63) 

Such scaling functions with zero moments are very useful for us from the point 
of view of time-frequency localization, because we have for Fourier component 



$(c<j) of them, that exists some C(ip) G R, such that for uj — > $(u;) = 1 + 
C(</?) | cj | 2r+2 (remember, that we consider r-regular multiresolution analysis). 
Using such type of scaling functions lead to super convergence properties for general 
Galerkin approximation [13]. Now we need some estimates in each scale for non- 
linear terms of type u \— > f(u) — f o u, where f is C°° (in previous and future parts 
we consider only truncated Taylor series action). Let us consider non regular space 
of approximation V of the form 

V = V q ® £ W s , (64) 

with Wj C Wj. We need efficient and precise estimate of / o u on V. Let us set for 
q G Z and u G ^ 2 (R) 

]J f q (u) = 2~i/ 2 J2 /(2 ?/2 < % <p q , k >) • <p q ,k (65) 

fcez 

We have the following important for us estimation (uniformly in q) for u,f(u) G 

#(" +1 ) [13]: 

\\P q (/(«)) -II /«(«) ||l 2 =0(2-(" + ^) (66) 

For non regular spaces (64) we set 

n/v(«)=n /»+ e ^n/ww (67) 

^= ? ,p-i 

Then we have the following estimate: 

ii^ (/(«)) -n/v(«)ii^=°( 2 " (n+1)ff ) ( 6g ) 

uniformly in q and V (64). 

This estimate depends on q, not p, i.e. on the scale of the coarse grid, not on 
the finest grid used in definition of V. We have for total error 

ll/(«) " II /v(«)H = H/(«) - P v(f( u ))h> + W P v(f( u ) - II /yH)IU 2 (69) 

and since the projection error in V: \\f(u) — Py (/(«)) \\l 2 is much smaller than 
the projection error in V q we have the improvement (68) of (66). In our concrete 
calculations and estimates it is very useful to consider our approximations in the 
particular case of c-structured space: 

P-i 
V = V q + J2 span{i/j jik , k G [2 {j ~ l) - c, 2 {j ' l) + c] mod 2 j } (70) 

j=q 



V VARIATIONAL WAVELET APPROACH 
FOR PERIODIC TRAJECTORIES 

We start with extension of our approach to the case of periodic trajectories. 
The equations of motion corresponding to Hamiltonians (from part II) may also 
be formulated as a particular case of the general system of ordinary differential 
equations d,Xi/dt = fi(xj,t), (i,j = 1, ...,n), < t < 1, where f i are not more than 
polynomial functions of dynamical variables Xj and have arbitrary dependence of 
time but with periodic boundary conditions. According to our variational approach 
we have the solution in the following form 

Xi (t) = Xi (0) + ]T A*V*W, Xi(0) = Xi(l), (71) 

k 

where Af are again the roots of reduced algebraical systems of equations with the 
same degree of nonlinearity and </?&(£) corresponds to useful type of wavelet bases 
(frames). It should be noted that coefficients of reduced algebraical system are the 
solutions of additional linear problem and also depend on particular type of wavelet 
construction and type of bases. 

This linear problem is our second reduced algebraical problem. We need to find 
in general situation objects 

oo 

AS.'.t=/n^;Wdx, (72) 

— oo 

but now in the case of periodic boundary conditions. Now we consider the procedure 
of their calculations in case of periodic boundary conditions in the base of periodic 
wavelet functions on the interval [0,1] and corresponding expansion (71) inside our 
variational approach. Periodization procedure gives us 

¥j,k( X ) = Yl Vhk{ x ~ I) ( 73 ) 

So, (p,ip are periodic functions on the interval [0,1]. Because ipj^ = ipj^' if k — 
/c'mod(2J), we may consider only < k < 2- ? and as consequence our multiresolution 

has the form M Vj — L 2 [0, 1] with Vj = span{(^j k } k -!r o [14]. Integration by parts 

i>o 
and periodicity gives useful relations between objects (72) in particular quadratic 
case (d = d\ + d 2 ): 

\di,d,2 / T\di \0,d,2+di k0,d *0,(i __ id (7A\ 

iV fcl,fc 2 — V l > 1V k 1: k 2 ' 7V fci,fc 2 — 1V 0,k 2 -ki — iV fc 2 -fci \'*J 

So, any 2-tuple can be represent by A^. Then our second additional linear problem 
is reduced to the eigenvalue problem for {A^} < A .< 2 i by creating a system of 2 Jr 



homogeneous relations in A. d and inhomogeneous equations. So, if we have dilation 
equation in the form (p(x) = V2J2k&z hk*-p(2x — k), then we have the following 
homogeneous relations 

N-1N-1 

A* = 2 d E E h m h e A d i+2k _ m , (75) 

m=0 e=o 



or in such form A\ d = 2 d \ d , where \ d = {A^} <fc<2J- Inhomogeneous equations 
are: 

E M d A d = d\2~ ]/ \ (76) 

where objects M d (\£\ < N — 2) can be computed by recursive procedure 

M d = 2 -i(2d+i)/2^d > Mk =< ^ ^ >= j2 ( k ] n k - j Mi, M l = 1. (77) 

So, we reduced our last problem to standard linear algebraical problem. Then 
we use the same methods as in part III C. As a result we obtained for closed 
trajectories of orbital dynamics described by Hamiltonians from part II the explicit 
time solution (71) in the base of periodized wavelets (73). 

VI VARIATIONAL APPROACH IN BIORTHOGONAL 

WAVELET BASES 

Now we consider further generalization of our variational wavelet approach. In 
[l]-[3] we consider different types of variational principles which give us weak solu- 
tions of our nonlinear problems. 

Before we consider the generalization of our wavelet variational approach to the 
symplectic invariant calculation of closed loops in Hamiltonian systems [3]. We also 
have the parametrization of our solution by some reduced algebraical problem but 
in contrast to the general case where the solution is parametrized by construction 
based on scalar refinement equation, in symplectic case we have parametrization 
of the solution by matrix problems - Quadratic Mirror Filters equations [3]. But 
because integrand of variational functionals is represented by bilinear form (scalar 
product) it seems more reasonable to consider wavelet constructions [15] which take 
into account all advantages of this structure. 

The action functional for loops in the phase space is [16] 



F( 7 )= fpdq- f H(t n (t))dt (78 

J*i Jo 



The critical points of F are those loops 7, which solve the Hamiltonian equations 
associated with the Hamiltonian H and hence are periodic orbits. By the way, all 



critical points of F are the saddle points of infinite Morse index, but surprisingly 
this approach is very effective. This will be demonstrated using several variational 
techniques starting from minimax due to Rabinowitz and ending with Floer ho- 
mology. So, (M, a;) is symplectic manifolds, H : M — > R, H is Hamiltonian, Xh is 
unique Hamiltonian vector field defined by 

w(X H (x),v) = -dH(x)(v), veT x M, x G M, (79) 

where u is the symplectic structure. A T-periodic solution x(t) of the Hamiltonian 
equations 

x = X H (x) on M (80) 

is a solution, satisfying the boundary conditions x(T) = x(0),T > 0. Let us 
consider the loop space Q = C°°(S 1 ,R 2n ), where S* 1 = R/Z, of smooth loops in 
R 2n . Let us define a function $ : Q — > R by setting 

a i ,1 

$(x) = / - < -Jx,x> dt - / H(x(t))dt, x en (81) 

jo 2 Jo 

The critical points of $ are the periodic solutions of x — Xh(x). Computing the 
derivative at x G Q in the direction of y G Q, we find 

&(x)(y) = -<S>{x + ey)\ e=0 = [ < -Jx - \jH{x),y > dt (82) 

ae Jo 

Consequently, $'(x)(y) = for all y G fi iff the loop x satisfies the equation 

-Jx(t) - ^H(x(t)) = 0, (83) 

i.e. x(t) is a solution of the Hamiltonian equations, which also satisfies x(0) = x(l), 
i.e. periodic of period 1. Periodic loops may be represented by their Fourier series: 

x(t) = E e k2nJt x k , x k eR 2 \ (84) 

fcez 

where J is quasicomplex structure. We give relations between quasicomplex struc- 
ture and wavelets in our second paper in this volume (see also [3]). But now we 
need to take into account underlying bilinear structure via wavelets. 

We started with two hierarchical sequences of approximations spaces [15]: 

. . . V-2 C V_i C V C Vy C V 2 . . . , . . . V_ 2 C VLi C V C Vy C V 2 . . . , (85) 

and as usually, Wq is complement to Vo in Vy, but now not necessarily orthogonal 
complement. New orthogonality conditions have now the following form: 

W ±V , W ±V Q , Vj±Wj, VjlWj (86) 



translates of ip span W , translates of ip span W . Biorthogonality conditions 
are 



<^jk,^fk' >= ip jk (x)il) fk >(x)dx = 5 kk >5 jy , (87) 

J ~ oo 

where ipjk(x) = 2?l 2 i\){2?x — k). Functions (p(x), <p(x — k) form dual pair: 

<ip(x-k),<p{x-£) >=6 kh <ip(x-k),4>(x-£)>=0 for Vfc, W. (88) 

Functions (p,<p generate a multiresolution analysis. tp(x — k), ip(x — k) are syn- 
thesis functions, <p(x — £), ip(x — £) are analysis functions. Synthesis functions are 
biorthogonal to analysis functions. Scaling spaces are orthogonal Jo dual wavelet 
spaces. Two multiresolutions are intertwining Vj + Wj = Vj+i, Vj + Wj = Vj + \. 
These are direct sums but not orthogonal sums. 

So, our representation for solution has now the form 

/(t) = E *;*(*), (89) 

where synthesis wavelets are used to synthesize the function. But bj k come from 
inner products with analysis wavelets. Biorthogonality yields 

b tm = f f(t)ij>e m (t)dt. (90) 

So, now we can introduce this more complicated construction into our variational 
approach. We have modification only on the level of computing coefficients of 
reduced nonlinear algebraical system. This new construction is more flexible. 
Biorthogonal point of view is more stable under the action of large class of op- 
erators while orthogonal (one scale for multiresolution) is fragile, all computations 
are much more simpler and we accelerate the rate of convergence. In all types 
of Hamiltonian calculation, which are based on some bilinear structures (symplec- 
tic or Poissonian structures, bilinear form of integrand in variational integral) this 
framework leads to greater success. 

VII VARIABLE COEFFICIENTS 

In the case when we have situation when our problem is described a system of 
nonlinear (polynomial)differential equations, we need to consider extension of our 
previous approach which can take into account any type of variable coefficients 
(periodic, regular or singular). We can produce such approach if we add in our 
construction additional refinement equation, which encoded all information about 
variable coefficients [17]. According to our variational approach we need to compute 
integrals of the form 



b tJ (t)^i) dl (2 m t - k 1 )( ( p 2 ) d %2 m t - k 2 )dx, (91) 

D 

where now bij(t) are arbitrary functions of time, where trial functions tpi, (p 2 satisfy 
a refinement equations: 

tpi(t) = J2 aik<Pi(2t - k) (92) 

feez 

If we consider all computations in the class of compactly supported wavelets then 
only a finite number of coefficients do not vanish. To approximate the non-constant 
coefficients, we need choose a different refinable function (p 3 along with some local 
approximation scheme 

(B e f)(x) := £ F t , k (f)<p 3 (?t - k), (93) 

qGZ 

where F^ k are suitable functionals supported in a small neighborhood of 2~ e k and 
then replace bij in (91) by Bebij(t). In particular case one can take a characteristic 
function and can thus approximate non-smooth coefficients locally. To guarantee 
sufficient accuracy of the resulting approximation to (91) it is important to have 
the flexibility of choosing tp s different from ipi,ip 2 . In the case when D is some 
domain, we can write 

K(t)\o= E b t ,(t) XD (2 e t-k), (94) 

0<fc<2< 

where xd is characteristic function of D. So, if we take <y? 4 = xd, which is again 
a refinable function, then the problem of computation of (91) is reduced to the 
problem of calculation of integral 

H(k 1 ,k 2 ,k 3 ,k 4 ) = H(k)= (95) 

/ ipA2H - k 1 )ip 3 (2 e t - k 2 )(pi 1 (2 r t - fc 3 V2 2 (2 s t - A; 4 )dx 

The key point is that these integrals also satisfy some sort of refinement equation: 
2-^H(k) = ]T b 2k . e H(e), ix = h + d 2 . (96) 

This equation can be interpreted as the problem of computing an eigenvector. 
Thus, we reduced the problem of extension of our method to the case of variable 
coefficients to the same standard algebraical problem as in the preceding sections. 
So, the general scheme is the same one and we have only one more additional linear 
algebraic problem by which we in the same way can parameterize the solutions of 
corresponding problem. 

Extended version and related results may be found in [l]-[6]. 
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